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Abstract 

Using Langevin dynamics simulations, we investi- 
gate the dynamics of chaperone-assisted transloca- 
tion of a flexible polymer through a nanopore. We 
find that increasing the binding energy £ between 
the chaperone and the chain and the chaperone 
concentration N c can greatly improve the translo- 
cation probability. Particularly, with increasing the 
chaperone concentration a maximum translocation 
probability is observed for weak binding. For a 
fixed chaperone concentration, the histogram of 
translocation time X has a transition from long- 
tailed distribution to Gaussian distribution with in- 
creasing £. t rapidly decreases and then almost 
saturates with increasing binding energy for short 
chain, however, it has a minimum for longer chains 
at lower chaperone concentration. We also show 
that T has a minimum as a function of the chaper- 
one concentration. For different £, a nonuniversal 
dependence of x on the chain length N is also ob- 
served. These results can be interpreted by char- 
acteristic entropic effects for flexible polymers in- 
duced by either crowding effect from high chap- 
erone concentration or the intersegmental binding 
for the high binding energy. 

Introduction 

The transport of biopolymers through a nanopore 
embedded in a membrane has attracted wide at- 
tention because it closely connects with polymer 
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physics and is also related to many crucial pro- 
cesses in biology, examples including the pas- 
sage of DNA and RNA through nuclear pores, 
the translocation of proteins through the endoplas- 
mic reticulum, as well as the viral injection of 
DNA into a host.- In addition, the translocation 
processes have been suggested to have potentially 
revolutionary technological applications, such as 
rapid DNA or RNA sequencing,-"^ gene therapy 5 
and controlled drug delivery. 6 

Polymer translocation through a nanopore faces 
a large entropic barrier due to the loss of a great 
number of available conformations, thus driving 
forces are introduced. Two important driving 
forces for translocation, both in experimental se- 
tups and in vivo, are provided by an electric field 
across the membrane and binding proteins (so- 
called chaperones). Translocation driven by an 
electric field has recently been investigated ex- 
tensively.-""— In this study, we focus on the lat- 
ter driving mechanism, which is responsible for 
the translocation of proteins-^- 5 - as well as the 
DNA translocation through membranes.— 1 ^ Up to 
now, there are also several theoretical and numer- 
ical studies^""— specifically devoted to the chain 
translocation in the presence of chaperones. 

In the pioneering study by Simon et al.,— the 
role of chaperones has been recognized as a Brow- 
nian ratchet where the effect of chaperone's bind- 
ing to the protein is to prohibit its backward dif- 
fusion through the pore and consequently speed 
up the translocation. This mechanism rectifies the 
dynamics of the protein compared with pure diffu- 
sion case. Later, Zandi et al.— investigated the 
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chaperone-assisted translocation of a stiff poly- 
mer using Brownian dynamics and proposed a new 
mechanism. At the moment when binding event 
occurs, the chain experiences a net force along its 
length, pulling the chain. This result demonstrates 
that the role of chaperone is quite different from 
that of a Brownian ratchet. In a more detailed 
theoretical investigation based on master equa- 
tion, Ambjornsson et ah— =2 examined chaperone- 
driven translocation of a stiff polymer. They iden- 
tified three limiting dynamical regimes accord- 
ing to binding situations: slow binding where ei- 
ther the chaperone concentration is low or binding 
strength is small (diffusive regime), fast binding 
but slow unbinding in which the chain cannot slide 
backwards (the irreversible binding regime), and 
fast binding and unbinding (the reversible binding 
regime). 

However, the aforementioned studies did not 
take into account the chain flexibility due to 
the difficulties in theoretical treatment. Particu- 
larly, the chaperone's intersegment binding will 
lead to the bending of the polymer, thus the is- 
sue of chain flexibility is clearly important for 
chaperone-assisted translocation. The basic ques- 
tions associated with this process are the follow- 
ing: (a) What's the effect of the chaperone concen- 
tration and the binding energy on the translocation 
probability? (b) How do the chaperone concentra- 
tion and the binding energy affect the translocation 
time? Is there an optimum chaperone concentra- 
tion or an optimum binding energy for transloca- 
tion? These questions are very complicated and 
still not clear. To this end, using Langevin dynam- 
ics we investigate chaperone-assisted translocation 
of a flexible polymer through a nanopore. 

Model and methods 

In our simulations, the polymer chain is modeled 
as a bead- spring chain of Lennard- Jones (LJ) par- 
ticles with the Finite Extension Nonlinear Elastic 
(FENE) potential. Excluded volume interaction 
between monomers is modeled by a short range re- 
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pulsive LJ potential: U]j(r) = 4£o[( 
e for r < 2 1 / 6 <7 and for r > 2 1 / 6 a. Here, a 
is the diameter of a bead, and £o is the depth of 
the potential. The connectivity between neighbor- 
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Figure 1 : (Color online) Schematic illustration of 
polymer translocation through a nanopore in the 
presence of chaperones depicted by blue particles. 
The width of the pore is w = 1.6a. 



ing monomers is modeled as a FENE spring with 
UFENE(r) — — 5&i?Qln(l — r 2 /Rq), where r is the 
distance between consecutive monomers, k is the 
spring constant and Rq is the maximum allowed 
separation between connected monomers. 
We consider a two-dimensional geometry as de- 

where the chaperones of 
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diameter a are modeled as mobile beads moving 
within the rectangular box with a pore of width 
w = 1.6 a formed by stationary wall particles of 
the same size. The wall particle interacts with 
the monomer and the chaperone by the repulsive 
LJ potential as shown above, while chaperone- 
monomer interaction is modeled by an attractive 
LJ potential with a cutoff of 2.5a and interaction 
energy £. In the Langevin dynamics simulations, 
each monomer is subjected to conservative, fric- 
tional and random forces, respectively. Namely, 
mfi = -V(U U + U F ene) - £, V; + Ff , where m is 
the bead's mass, ^ is the friction coefficient, v ; is 
the bead's velocity, and Ff is the random force 
which satisfies the fluctuation-dissipation theo- 
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In our model, the LJ parameters £o, a and bead 
mass m fix the system energy, length and mass 
scales, leading to the corresponding time scale 
tu = (mo 2 /£o) 1 / 2 . Each bead corresponds to a 
Kuhn length of a polymer, so we choose a ~ 1.5 
nm and the bead mass m fa 936 amu.— We set 
ksT = 1.2£o, which means the interaction energy 
£o to be 3.39 x 10~ 21 J at actual temperature 295 
K. This leads to a time scale of 32.1 ps. - 3 The 
dimensionless parameters in our simulations are 
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chosen to be R = 1.5, k = 30, = 0.7. Then 
the Langevin equation is integrated in time by the 
method proposed by Ermak and Buckholz.— Ini- 
tially, the first monomer is placed at the trans side 
with one unit length to the pore center. Then, 
the remaining monomers and chaperones are un- 
der thermal collisions described by the Langevin 
thermostat to reach equilibrium state of the sys- 
tem. Typically, each simulation data is the aver- 
age result of 1000 successful translocation events 
to minimize statistical errors. 

The current simulation method was also used 
to investigate polymer translocation driven in the 
presence of a cross-membrane electrical poten- 
tialil^ 

Results and discussions 

Based on master equation, Ambjornsson et al.—^ 
investigated the dynamics of chaperone-driven 
translocation of a stiff polymer. Although poly- 
meric degrees of freedom of the translocating 
chain is neglected and the system is assumed to 
be close to equilibrium, their results still can shed 
light on the translocation of the flexible polymer. 

For the chaperone binding to the binding site, the 
probability that a binding site is occupied depends 
on the chaperone concentration cq and the bind- 
ing energy £. By calculating the binding partition 
function, they obtain a dimensionless parameter 
K = coK eq as a relevant measure of the effective 
binding strength,^22, where the equilibrium bind- 
ing constant K eq = VQexpiejksT) with vo being 
the typical chaperone volume. For univalent bind- 
ing, the equilibrium probability that a binding site 
is occupied is 

PZ = K/{l + K). (1) 

For stiff polymers, the force acting on polymer, 
Fbind (s) , in unit of kgT /a for reversible binding of 
chaperones to the translocating polymer is— 

F bind (s)=ln{l + K). (2) 

where s is the translocation coordinates. This re- 
sults indicates that the force increases with the 
chaperone concentration and the binding energy 
and is independent of s. 



For flexible polymers, we don't expect that the 
effective binding strength fc and the probability 
that a binding site is occupied Pl% are the same 
as those for stiff polymers. Due to the chain flexi- 
bility, one chaperone can attach to several binding 
sites (intersegmental binding) for the high chaper- 
one concentration and binding energy as observed 
in the simulation, leading to chain folding. But k 
and Pqcc should also increase with the chaperone 
concentration and the binding energy. 

In addition, for flexible polymers entropic ef- 
fects from cis and trans sides give an additional 
contribution to the total force. As noted by Amb- 
jornsson et al.,~ the binding force and the entropic 
force are additive if the chaperone binding is inde- 
pendent of the curvature of the polymer. For flex- 
ible polymers, we can write the total force acting 
on the translocating polymer as 

F(s) = F bind (s) -F trans , e (s,c ,e) -F cis , e {s). (3) 

Here, F b i nc i(s) is the force from the binding which 
should increase with the chaperone concentration 
and the binding energy, although its explicit ex- 
pression may be different from ?? for high bind- 
ing energy. F trans ^ e (s,co,£) is the entropic force 
from the trans side, which also increases with 
the chaperone concentration and the binding en- 
ergy, resulting from either crowding effect induced 
by high chaperone concentration or the interseg- 
mental binding. The entropic force from the cis 
side F c j S; e(s) plays an important role only at the 
beginning of the translocation, and it is negli- 
gible compared with F tranS)e (s 1 co 1 e) particularly 
for high chaperone concentration and binding en- 
ergy. Therefore, the translocation dynamics is 
dominated by the interplay between Ft, inc i(s) and 

Ftrans,e(s i CQi C) • 

In our simulation, the chaperone concentration 
is proportional to the ratio of the number of the 
chaperone (7Y C ) to the constant area of the rectan- 
gular box (64 x 64). In the following, we use N c to 
stand for the chaperone concentration. 

The effect of e and N c on the transloca- 
tion probability 

The translocation probability, P trans , refers to the 
ratio of the successful translocation events to the 
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Figure 2: The translocation probability as a func- 
tion of (a) the binding energy £ for different chap- 
erone concentration (~ N c ) and (b) the chaperone 
concentration for different binding energy. The 
chain length N = 64. The inset shows situations 
for weak binding. 



whole attempts at given physical parameters in 
the simulation. With increasing £, P tra ns increases 
rapidly first, and then slowly approaches satura- 
tion at larger £ for higher chaperone concentra- 
tion, while it continuously increases for lower N c , 
as shown in |[figure][2][]2[ a). 

The observed P tra ns can be well understood by 
taking into account the force in ?? as a function of 
N c and £. The force exerted on the chain by chap- 
erone's binding, Fbi„d(s), increases with increas- 
ing the chaperone concentration and the binding 
energy £. For high binding energy £, F bind {s) is the 
dominant term in ?? for the translocation process. 
Meanwhile, unbinding is slower, as a consequence 
the chaperone acts as a ratchet which would ef- 
fectively prohibit the chain's backward motion out 
of the pore. Qualitative, P tra ns shows the similar 
behaviors as the probability that a binding sites is 
occupied, see 
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[figure] [2][]2f b) shows that P tr ans also goes up 
rapidly first with increasing Af c and then slowly 
approaches the saturation at higher N c for strong 
binding. The reason for the increase of Ptrans 
with N c is that the probability of chaperone's 
binding to the chain rises with increasing _/V c . 
However, for weak binding P tra ns shows a max- 
imum with increasing iV c , as shown in the inset 
of [figure][2][]2| ^b). This unexpected decrease for 
high chaperone concentration stems from crowd- 
ing effect at the trans side. Obviously, Fb ina -(s) 
increases with A^ c , however, Ft ra ns , e (5,co,£) also 
increases with growing Af c due to fact that the 
chain has to overcome greater entropic barrier for 
translocation and becomes dominant. Moreover, 
for higher chaperone concentration, the collision 
frequency between chaperones increases, which 
leads to higher unbinding events particularly for 
weak binding. The interplay of these factors re- 
sults in the maximum of Ptrans- 



Distribution of the translocation time 

We also checked that the binding energy £ has an 
obvious effect on the shape of the histogram of the 
translocation time, as shown in |[figure] [3][]3| For 
iV = 64 and A^ c = 650, the distribution of translo- 
cation time is asymmetric with a long exponen- 
tial tail for weak binding £ = 1.3, while for strong 
binding £ = 5.0, it nearly approaches a Gaussian 
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Figure 3: (Color online) The distribution of 
translocation time for different binding energy un- 
der N c = 650. The chain length is N = 64. 



distribution. The reason is that with increasing £ 
the force from binding F^^s) greatly increases 
and is the dominant term in the total force F(s). 
For driven translocation by cross -membrane elec- 
tric field, we also observed similar distributions 
for very weak and strong driving forces, respec- 
tively.— 

Translocation time as a function of e 

and N c 



[figure][4][]4^a) shows T as a function of £ for 
N = 64 and different N c . We find that x decreases 
with increasing £ (at least for £ < 8.5). It is not 
difficult to understand the overall fall of x with £ 
from the perspective of the magnitude of F(s) and 
the ratchet mechanism that prevents back-sliding 
motion as noted above. Increasing £ leads to the 
increase of Ft>i n d(s), and is beneficial to the for- 
ward motion of the chain, speeding up the translo- 
cation. What draws our attention most is the no- 
table discrepancy in the decay rate of x, i.e., x 
decreases quite rapidly with £ for £ lower than 
a critical binding energy e c , but much slowlier at 
£ > £ c regimes. Specifically, for £ < e c , the bind- 
ing frequency is higher but the force F bind {s) from 
the binding is smaller compared with the case of 
strong binding. The latter dominates the translo- 
cation dynamics, leading to the observed behav- 
ior. For £ > £ c , translocation is slowed due to 




Figure 4: (Color online) Translocation time as a 
function of the binding energy for (a) N = 64 and 
different chaperone concentration, and (b) N c = 30 
and different chain length N. 
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the competition between Ft, inc i(s) and the entropic 
force F trans>e (s,co,e), leading to slower increase of 
F(s) with increasing £. On the one hand, Fj ?iru i(s) 
increases with £. On the other hand, intersegmen- 
tal binding between translocated monomers results 
in larger entropic force F trans ^ e {s,CQ,e) due to the 
loss of chain conformations. 

However, for longer chain length N > 128 and 
lower chaperone concentration N c = 30 as shown 
in [figure] [4] []4| fo), we observed an optimum £ for 
translocation, namely X has a minimum as a func- 
tion of £. The time t^iff for a polymer of length 
N to diffuse a distance of the order of the binding 
site length o is x^iff ~ o 2 / D ^ N^o 2 / {k B T) with 
the diffusion constant of the chain D = ksT / (NE, ) . 
The average distance between chaperones in solu- 
tion is R c = Lj \fW c , where L — 64 is the length 
of the simulation box. It is suffice for a chap- 
erone to diffuse a distance of the order of R c for 
any one chaperone to attach the binding sites (pro- 
vided the binding energy is sufficiently high), and 
this time is x unocc = R 2 /4D C = $e 2 L 2 /(k B TN c ) = 
34. \^o 2 /(ksT), where the diffusion constant of 
chaperones D c = kgT /E, due to the same size for 
a chaperone and a monomer. Due to diffusion 
through a nanopore, x^i// has a large prefactor.— 
Thus, Xdiff is much longer than x unocc for larger N, 
which indicates that it is possible for chaperones to 
bind the chain very soon by taking into account the 
range of interaction of the cutoff 2.5a used in the 
simulation. For chain length ./V < 128, the chap- 
erone can almost cover all the binding sites for 
large £ by the intersegmental binding, thus there 
is no minimum for translocation time. But for 
longer chains, all chaperones have become bound 
and there are no free chaperones left before com- 
pleting the translocation, leading to the increase of 
the translocation time. 

To understand this unexpected behavior, we have 
also investigated the distribution for waiting (resi- 
dence) time of bead s, which is defined as the time 
between the events that the bead s and the bead 



s+1 exit the pore, [figure] [5] []5 shows the waiting 
time distribution for N = 256, iV c = 30 and differ- 
ent binding strength £ = 2.5 and 7.5, respectively. 
For £ = 2.5, the waiting time increases rapidly first 
and then approaches the saturation very soon, in- 
dicating the feature of fast binding and unbinding. 
However, for £ = 7.5 the waiting time always in- 
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Figure 5: Waiting time distribution for N = 256, 
N c = 30 and different binding strength £. 

creases except for last several monomers. What's 
more, the translocation is faster for £ = 7.5 than 
for £ = 2.5 when s < 125, however, it is much 
slower for £ = 7.5 than for £ = 2.5 when s > 125. 
All together, the translocation time for £ = 7.5 
is much longer than that for £ = 2.5. The rea- 
son for observed behavior for s > 125 is that, for 
£ = 7.5 > e c and low chaperone concentration, al- 
most all chaperones keep bound to the front part 
of the chain, rendering that there is no free chap- 
erones for latter new-emerging segments and the 
corresponding translocation becomes diffusive.— 




Figure 6: (Color online) The influence of chap- 
erone concentration on the translocation time for 
different £. The chain length is N = 64. 



[figure] [6] [] 6 shows x as a function of iV c for 



different £. We find that x initially decreases 
and subsequently goes up with increasing N c . As 
stated above, growing /V c leads to an increase in 
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Fbind(s) as well as F tranS;e (s,c ,e). Considering 
that Ffri nc i(s) is the dominant term in ?? for initial 
increasing of N c , the total force F(s) increases, re- 
sulting in the decrease of x. Besides, the proba- 
bility that the binding site just passing through the 
pore exit also increases with growing N c , which is 
favorable to the formation of the ratchet. 

However, overfull chaperones gives rise to 
the crowding effect and the entropic force 
F tra ns,e( s i c 0i £ ) prevails in the competition with 
Fbi?id( s ), hindering the translocation and result- 
ing in the increase in x instead of continuous fall. 
Moreover, the corresponding chaperone concen- 
tration of the minimum shifts to higher N c value 
with increasing £. This is due to the fact that with 
increasing £, the superiority of F tranSje (s,CQ,£) 
works for higher iV c and the formed ratchet could 
prevent the back-sliding motion of the chain more 
effectively. Needless to say, it is conducive for the 
translocation of the chain through the nanopore 
and thus postpones the advent of the minimum 
point. 

Translocation time as a function of the 
chain length 



length N is an important measure of the underly- 
ing dynamics, x ~ N a with a being the scaling 
exponent. Our previous two dimensional electric- 
field driven translocation simulations- 3 *^ show 
there was a crossover from a = 1V2D for fast 
translocation to a = 1 ± V2d for slow transloca- 
tion, where V2d = 0-75 is the Flory exponent for 
a self-avoiding chain in two dimensions. 34 Most 
recently, we further find that for fast translocation 
processes a = 1.37 in three dimensions, while it 
crosses over to a = 1 + V30 with V30 = 0.588 for 
slow translocation, corresponding to weak driving 
forces and/or high friction.— 



[figure] [7] []7 shows the x as a function of N for 
moderate chaperone concentration N c = 650 and 
different £. Obviously, the scaling exponent a de- 
pends significantly on e : it initially decreases from 
1.97 ±0.07 to a minimum 1.41 ±0.01, following 
by a slight increase with increasing £. For higher 
£, the force from the binding F bind (s) is greater 
and dominates the translocation dynamics, leading 
to a ~ 2V2D for fast translocation process as the 
electric-field driven translocation.- 3 -^ 3 - 3 - Decreasing 
£ to £ = 1.3, a increases to 1.97. If decreas- 
ing £ further, it may access unbiased translocation 
regime with a = N l+2v2D .— 
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Figure 7: (Color online) Translocation time as a 
function of the chain length N for different bind- 
ing energy £ under fixed chaperone concentration 
iV c = 650. 

Similar with the electric-field driven case, the 
presence of chaperones only in the trans side 
could also induce a chemical potential difference 
between the two sides of the membrane. The 
scaling of the translocation time x with the chain 



Conclusions 

Using Langevin dynamics simulations, we inves- 
tigate the dynamics of chaperone-assisted translo- 
cation of a flexible polymer through a nanopore. 
We find that increasing the binding energy £ be- 
tween the chaperone and the chain and the chaper- 
one concentration can greatly improve the translo- 
cation probability. Particularly, with increasing the 
chaperone concentration a maximum translocation 
probability is observed for weak binding. For a 
fixed chaperone concentration, the histogram of 
translocation time x has a transition from long- 
tailed distribution to Gaussian distribution with in- 
creasing £. x rapidly decreases and then almost 
saturates with increasing binding energy for short 
chain, however it has a minimum for longer chains 
at lower chaperone concentration. We also show 
that x has a minimum as a function of the chap- 
erone concentration. For different £, a nonuni- 
versal dependence of x on the chain length Af is 
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also observed. These results can be interpreted 
by characteristic entropic effects for flexible poly- 
mers induced by either the crowding effect from 
high chaperone concentration or the intersegmen- 
tal binding for the high binding energy. 

Generally, chaperones have a size larger than 
that of a monomer, giving rise to the "parking lot 
effect" as observed by previous studies,—"— lead- 
ing to a less efficient translocation: after binding 
of a chaperone to the chain close to the pore exit, 
the chain firstly needs to diffuse by a chaperone 
size distance, before next binding event occurring. 
In addition, we have also assumed that the binding 
energy is the same along the chain. However, pro- 
teins and nucleic acids consist of heterogeneous 
sequence of aminoacids, bases or base pairs, re- 
spectively. It has been found that the chain het- 
erogeneity is important in translocation dynamics 
for stiff polymers. 27 In the future studies, it would 
be interesting to investigate the effects of the size 
difference, the changes in chain flexibility and the 
chain heterogeneity along the chain on the translo- 
cation dynamics. 
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